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, Abstract 

^vq ■ We study the diffractive production of qqg jets with large transverse momenta in the region 

C*~) \ of large diffractive masses (small j3). Cross sections for transverse and longitudinal photons are 

obtained in the leading log l/xp and log 1/(3 approximation, keeping all powers in log kf /Q 2 . 
| We perform a numerical study and illustrate the angular distribution of the three jets. We 

also estimate the integrated diffractive three jet cross section and compare with the dijet cross 
section obtained before. 
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1 Introduction 



The diffractive dissociation in deep inelastic scattering lies on the interface between perturbative and 
non-perturbative QCD, and one expects the diffractive cross section to contain both hard and soft 
contributions. As a hrst experimental hint in this direction, one might take the observation that the 
energy dependence of the diffractive cross section, i.e. the measured intercept of the Pomeron flux 
factor tends to be larger than that of the pure soft Pomeron Jl], [|]. As a possible interpretation of 
this observation, the contribution of the "perturbative Pomeron", the gluon structure function, may 
turn out to be substantial: the observed rise of the gluon structure function at small - x enhances 
the hard part of the diffractive final state and thus leads to an effective Pomeron intercept which 
is bigger than that of the soft Pomeron. 

In a first attempt to provide a description of the diffractive cross section which takes into account 
both hard and soft contributions a parameterization has been suggested || and successfully fitted 
to both ZEUS and HI g g data. As a main result, it was shown that an acceptable description 
of the data requires (at least) three pieces: the diffractive production of transverse qq and qqg 
states and the diffractive production of longitudinal vector mesons. The fits also confirmed that 
the effective intercept of the Pomeron tends to be larger than the soft Pomeron. 

In order to obtain further insight into DIS diffraction and, in particular, to understand the origin 
of the higher Pomeron intercept, one has to open the diffractive final state and to isolate those parts 
which can be attributed to the hard Pomeron. Apart from the longitudinal vector particles which 

1 Supported by the TMR Network "QCD and Deep Structure of Elementary Particles" 
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have been suggested already several years ago [p|, diffractive final states consisting of hard jets only, 
are the most promising candidates. As a first step, the production of quark- antiquark jets (Fig. |l|) 
with large transverse momenta has been discussed by several groups ||, ^, 0. To summarize the 
main result, the energy dependence is given by the square of the gluon structure function: 

da ~ [-L XF g(x F , 4/(1 - I3))f (1-1) 

where kx denotes the transverse momentum of the quarks in the photon - Pomeron rest frame, and 
(5 = Q 2 / (Q 2 + M 2 ) is related to the mass M of the diffractive system. From eq. (1.1) one concludes 
that for small masses ((3 near 1) the cross section will be dominated by this "hard Pomeron", either 
in the production of vector particles or of hard gg-jets. For larger diffractive masses the contribution 
of hard qq jets is smaller but not negligible. Another interesting feature is the azimuthal dependence 
of the qq cross section. It was found || that, in the photon - Pomeron rest frame, the two jets prefer 
to lie perpendicular to the scattering plane (defined by the photon - Pomeron direction and the 
outgoing electron momentum). This has to be compared with the photon - gluon fusion mechanism 
in the usual DIS final state, where the two outgoing quarks prefer to lie inside the scattering plane. 
The experimental verification of these theoretical predictions turns out to be rather difficult: to 
make the calculation applicable, the jet final state is not allowed to contain any soft Pomeron 
remnants; the presence of such pieces changes the situation rather radically. To exclude such a soft 
remnant requires high statistics of the data and a rather good resolution in the jet algorithm. First 



results have been reported in [10] . 
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Figure 1: Dijet Production. The exchanged gluons couple to the quarks in all possible 
ways. 

For not so small diffractive masses (/? away from 1) one expects the production of extra gluons 
to become essential. The simplest final state consists of the qq pair and one gluon; in the region of 
large diffractive masses (small (3) one expects a separation in rapidity between the qq system and the 
gluon (Fig. ID) (not to be confused with the rapidity gap between the diffractive qqg system and the 
outgoing proton). In order to isolate the "hard Pomeron" contributions one, again, wants to study 
final states which consist only of jets with large transverse momenta (and no further soft remnants). 
In this paper we present the calculation of this three-jet cross section, restricting ourselves to the 
region of large diffractive masses (small j3). We shall work in the leading - log(l//3), log(l/xip) 
approximation, retaining all powers in k^/Q 2 . This approximation includes also higher-twist terms 
of the diffractive cross section. Diffractive production of qqg final states in another kinematic region 



has been studied before [11, 12 1: the gluon transverse momentum has to be smaller than those of 



the quarks. The results of [12| extend into the large /3-region. 
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Figure 2: Three jet Production. The exchanged gluons couple to the quarks in all possible 
ways. 

There are several questions to which our jet analysis gives an answer. Apart from a quantitative 
estimate of the jet rates (as a function of the imposed kinematical cuts) our analytic formula allows 
a first step towards an analysis of the event shape of diffractive final states (e.g. thrust). In 
particular, we find that in the region of medium small (3 (where our formula is applicable) the 
qqg final state provides a substantial contribution to the two-jet cross section. Our framework 
also allows to address the transition from hard jets to the aligned gluon configuration: in a future 
step one might model a two gluon amplitude which allows to extrapolate our formula into the 
soft Pomeron region. Such a model (an explicit parameterization was suggested in ref. ||13||) would 
allow to describe both the hard and the soft region of diffractive qqg production. However, a severe 
limitation is our restriction to the small-/? region. In order to cover the full (3 range we need a 
complete LO calculation of three parton final states in DIS diffraction which includes a NLO two 
jet calculation. 

This paper will be organized as follows. In section 2 we describe our calculation and present 
the result for the three jet production cross section. Section 3 contains the transformation into the 
impact parameter space and in section 4 we discuss a few particular kinematic regions for which 
our formula allows analytic approximations. In section 5 we perform a first numerical analysis, 
calculate integrated cross sections and show a few plots. In the final section we summarize our 
results and make a few concluding remarks. 

2 The Differential Jet Production Cross Section 

Our notations are illustrated in Fig. ^. As usual, q and p denote the four momenta of the photon 
and the proton, respectively, t = r 2 is the square of the momentum transferred from the proton 
to the diffractive system. The total energy is denoted by s and the sub-energy of the 7*-proton 
system with W 2 . The scaling variables are x = Q 2 /2pq, xjp = (Q 2 + M 2 )/(Q 2 + W 2 ) , and 
P = Q 2 /(Q 2 + M 2 ). For simplicity we restrict ourselves to the forward direction, t = 0, such that 
the four momentum r = xpp has no transverse component. We use the Sudakov parameterization 

ki = and + (3ip + kiT, k 2 = -k 2 T , (2.1) 

where q' = q + xp. We work in a reference frame where photon and proton momenta are collinear 
(along the z - direction); the scattering plane (the x — z plane) is defined by the incoming and 
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Figure 3: Notations for the diffractive production of a qqg system. 



outgoing electron. Azimuthal angles refer to this scattering plane. In order to justify the use of 
perturbation theory we start in the "hard region": all three final state partons have large transverse 
momenta. Later on we will show that the more precise conditions imply that the virtualities of the 
exchanged quark and gluon have to be large. This leads to: 

Qo < k 2 

Ql < — { ai (l- ai )Q 2 + k 2 } 

Ql < -^Ml-ai)Q 2 + (ki+k 2 ) 2 ] (2.2) 

1 — Oil 

Our calculation will be done in the leading - log W 2 , leading - log M 2 approximation and applies 
to the region Q 2 <C M 2 <C W 2 . As usual, in this approximation we have a2 <C ai, 1 - ai and 02 
close to xjp. For large values of logQ 2 /k 2 our approximation reduces to the double logarithmic 
approximation (DLA). The invariant mass M of the diffractive system can be expressed (in this 
leading log approximation) as 

M 2 = W 2 (xjp-x) . (2.3) 
The mass of the qq subsystem is given by m 2 ^ = m 2 — k\ with 

m 2 = M 2 — [5 2 W\ 02 = ■ (2.4) 

For Q 2 <C M 2 (i.e. small (3) we expect the invariant mass of the qq pair to be much smaller than 
M 2 (of the order of Q 2 ). From the mass shell condition of the upper quark line (with momentum 
r + q — k\ — k 2 ) we have the relation: 

= a 1 (l-a l -a 2 )(M 2 - l3 2 W 2 )-{l-a 1 -a 2 )'k 2 -a 1 (k 1 + k 2 ) 2 

w «i(l - ai)m 2 - (1 - Qi)k 2 - «i(ki +k 2 ) 2 (2.5) 

We begin with the calculation of the j*p cross section. The diagrams we have to calculate are 
shown in Fig. ||. Since we consider the square of the scattering matrix element and, moreover, take 
the discontinuity in W 2 , all quark lines between the leftmost and the rightmost gluon are on-shell; 
the same holds for the produced gluon below the quark loop. This leads to ^-functions which will be 
used to integrate over all Oj and fy, except for (5 2 . We begin with the transverse cross section. Let 
e x = (0, 1, 0, 0) = (0, e x , 0) and e y = (0, 0, 1, 0) = (0, e y , 0) be the transverse polarization vectors of 
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the photon (it should be noted that e x lies in the scattering plane). Contracting the Lorentz indices 
of the two external photons with the transverse polarization vectors e\ and e 2 (each of them can 
be either e x or e y ), the sum of the diagrams in Fig. || leads to a contribution of the form 



2W A 5 [{r + q-h- k 2 f) [4a? (1 - ai) 2 e H e 2i M a Mj, - ai(l - a x )ei • e 2 MyM a J (2.6) 

(here / = 1,2 refer to the transverse components; repeated indices are summed over, and the 
prime at the second factor indicates that it belongs to the two-gluon amplitude with momentum 
1', attached to the lower two gluons on the rhs of the diagrams in Fig. 4). We have exhibited the 5 
function of the upper quark line which leads to the condition eq. (2.4) and will be used to express 
«i in terms of M, /3 2 , ki, and k 2 . The expression for M™ will be discussed below. Contracting this 
with the lepton tensor -^L^e^e 1 ^ and summing over the two possible transverse polarizations we 
obtain 



2W A 5 ((r + q-h- k 2 f) 



1 + (1 o v)2 A ai {\ - a x )[a\ + (1 - a^MaM'a - 4(1 - y)a\(l - a 1 f(M xl M[ l - M 2l M' 2l ) 



(2.7) 



In the second line one recognizes two terms. The first one (proportional to My AT',) corresponds to 
the sum over the two transverse polarizations, the second one (proportional to MijMy — M^M^, 
Mi; = e x iMn, M 2 i = e y iMn) to the difference. As long as we do not integrate over the azimuthal 
angle of the outgoing electron, the two transverse polarizations of the photon come with different 
weights, and that is why we write our result in terms of both the sum and the difference of the two 
polarizations. 

For the longitudinal cross section we arrive at the following structure: 



2W A 5 [(r + q-h- k 2 ) 2 ) a\(l - a 1 fQ 2 M l M[. (2.8) 

Together with the lepton tensor we obtain 

2W A 5 [(r + q-h- k 2 ) 2 ) 2(1 - y)a\(l - a 1 fQ 2 M l M[. (2.9) 

Finally, we have the interference between the transverse and the longitudinal photon. The 
expressions analogous to eq. (2.5) and eq. (2.6) are given by 



2W A 5 Ur + q-h- k 2 ) 2 j 2^Q 2 a\{l - ai) 2 (l - 2a 1 )M n M' l . (2.10) 

and 

2W 4 5 [(r + q-h- k 2 ) 2 ) yi^y(2 - y)^Q 2 al(l - ai) 2 (l - 2a x )M ll M[. (2.11) 

Before we turn to the calculation of the My and Mi, we write down the final formula for the 
process e + p — * e' + qqg + p. The diffractive final state is characterized by the invariant mass M, 
the longitudinal momentum fraction (3 2 of the outgoing gluon, and by the transverse momenta ki, 
k 2 of the outgoing quark and gluon, respectively (in our notation, k 2 = — k 2 -. > 0). Starting from 
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l + r I I' I'+r 
a 




c 

Figure 4: Diagrams for the cross section of diffractive qqg production. We only show the 
diffractive system: at the lower end the two-gluon systems on both sides couple to the 
gluon structure function of the proton. At the upper end, the gluons couple to the quark 
loop in all possible ways. We do not show the vertical discontinuity line which divides 
each diagram into left and right hand pieces. 
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eq. (2.5) - eq. (2.10) and evaluating the final state phase space integrals we obtain: 



da e D p 



a, 



dydQ 2 dM 2 dm 2 d 2 k 1 d 2 k 2 dt [t=0 
l + (l-y) 2 



daj) P T . 



dM 2 dm 2 d 2 k 1 d 2 k 2 dt 



|t=o 



yQ 2 ir 

2(1 - y) 



daj) P T _ 



+(i-y)- 



d<rV L 



dM 2 dm 2 dk l d 2 k 2 dt\ t=0 
where the ^y*p cross sections are defined as: 



+ (2-y)y/T=y 



dM 2 dm 2 d 2 k 1 d 2 k 2 dt {t=0 
daft 



dM 2 dm 2 d 2 k 1 d 2 k 2 dt 



t=o 



(2.12) 



a<J D,T+ 


9 1 




dM 2 dm 2 d 2 k 1 d 2 k 2 dt\ t=0 


128tt VS(M 2 - m T 
■MuM-i 


m 2 


daj) P T _ 


9 1 




dM 2 dm 2 d 2 k 1 d 2 k 2 dt lt=0 


128tt \f~S(M 2 - m T 


m 2 




{MuM'u - 
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dM 2 dm 2 d 2 k 1 d 2 k 2 dt lt=0 


128tt y/S{M 2 - m 2 ] 


m 2 
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dM 2 dm 2 d 2 k 1 d 2 k 2 dt ]t=Q 


128tt V5(M 2 - m T 


m 2 



E9 3 
e f a em a s 



a\ + (1 — ai) 2 ai(l — «i) 



f 



(2.13) 



(2.14) 



53 eja em a^4a? (1 - «i) 3 Q 2 M z M/ (2.15) 
/ 



53 e)a em ala\{l - «i) 2 (l - 2a x ) 
/ 



/Q 2 [M u M/ + MiM'u] . 



Here 



1 + 



k\ (ki + k 2 ) 



k 2 
K i 



(2.16) 



(2.17) 



and the prefactors containing ct\ have to be re-expressed using the ^-function condition eq. (2.4). 
As we have discussed before, the transverse cross section has two terms. The first one corresponds 
to the sum over the two transverse polarizations, the second one to the difference. This second term 
of the transverse cross section, as well as the interference cross section, introduce the azimuthal 
dependence of our cross section. 

It is instructive to consider the kinematics of the three partons in the final state. Within our 
approximation, we have to restrict ourselves to the forward direction of the Pomeron: in the 7*- 
1P CM system the gluon is emitted in the direction of the Pomeron. The kinematic boundaries 
of the transverse momenta ki and k 2 are determined by the condition S > (eq. (2.17)). The 
allowed region is illustrated in Fig. [5], where we have plotted the two quark transverse momenta, 
normalized to the invariant mass square m 2 of the qq subsystem: k 2 /m 2 and (ki + k2) 2 /m 2 . The 
curve S = shows that the quark momenta are essentially restricted by m 2 . For the particular case 
of small gluon transverse momenta k| <C k 2 , we are on the diagonal line k 2 ~ (ki + k 2 ) 2 , and the 
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(h + k 2 ) 2 /m 2 




k 2 /m 2 



Figure 5: Kinematic Boundaries. The drawn line belongs to S = 0. 



physical region is restricted to k 2 < \m 2 . For a given set of transverse momenta, the longitudinal 
momentum fraction a.\ is given by one of the two values: 



'f 1+ ^-(y*>w 

2 \ m z , 



(2.18) 



(02, Pi, and /?2 follow from the mass shell conditions k\ = 0, k\ = 0). The limiting values ai = 
and 1 — ai = lie on the axis k 2 = and (ki + k 2 ) 2 = 0, respectively, i.e. they correspond to final 
states where one of the quark lines becomes soft. For small k| (i.e. on the line kf = (ki + k2) 2 ) 
eq. (2.17) simplifies to 



ai 




and 



(2.19) 



ai(l — «i) = k^/m 2 . 



For the angle A between the two quarks we find 



1 — cos A 



(k x + a!k 2 ) 2 /M 2 



[a 2 + k 2 /M 2 ][(l - ai y + (ki + k 2 ) 2 /M 2 ] 



(2.20) 



(2.21) 



In the denominator, the a\ terms will dominate, as long as we stay away from the values a>\ = 0, 1, 
the numerator is (at most) of the order m 2 /M 2 and the opening ang le A = 0(m 2 /M 2 ). Therefore, 
for small m 2 /M 2 the opening angle will be small (of the order A = 0(m 2 /M 2 )), and we expect, in 
the 7* — IP CM system, the gluon and the gg-pair to come out as two back-to-back jets. The qq - pair 
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produces a broader jet, where the opening angle is given by eq. (2.21) and will approximately be 
of the order of 0(m 2 /M 2 ). When m 2 gets larger, the qq - pair jet becomes broader and eventually 
can be resolved into two separate jets. The angle of the axis of the two-jet configuration is given 
by the direction of the outgoing gluon. If # 2 denotes the angle between the gluon and the proton 
momentum, we find: 

l-cosfl 2 = k 2 /M 2 (2 22) 

2 k 2 /M 2 + (l-m 2 /M 2 ) 2 ' 1 ' ' 

For m 2 <C M 2 the gluon transverse momentum k| will be much smaller than M 2 , the angle will 
be small, i.e. the outgoing gluon moves in the direction of the proton with an angle of the order 
0(k 2 /M 2 ). In the strong ordering region k 2 <C k 2 , eq. (2.21) can be approximated by 

1 — cos A m 2 1 

2 = M 2 [a x + (1 - ai)m 2 /M 2 ][(l - ai ) + ai m 2 /M 2 } ( ' ' 

In section 5 we present results of a numerical study of our cross section formula which illustrate 
the geometry of the three-jet configuration. 

Let us now turn to the calculation of the Mn and M/ . Each of them is of the form 

M U = J ^T( X]P ,l 2 )T u , (2.24) 

(and an analogous expression for Mi), where T denotes the unintegrated gluon structure function 

rQ 2 

/ dl 2 F(x,l 2 ) =xg(x,Q 2 ) (2.25) 
JQI 

and models the perturbative Pomeron. In the fermion loop at the top of the diagrams in Fig. || 
it is understood that the gluons couple to the quarks in all possible ways. Since we consider the 
square of the scattering matrix element and take the discontinuity in W 2 , all quark lines between 
the leftmost and the rightmost gluon are on - shell; this leads to 5 - functions which will be used 
to integrate over all a, and except for fo- We begin with the transverse cross section. The 
diagrams with two gluons lines attached to the fermion loop (Fig. ||a) lead to: 

(a) _ / k 1+ k 2 kg \ / k 2 + l k 2 \ 

il ~ U(k 1+ k 2 ) J D(k 1 )JA(k 2 + l) 2 k 2 J ; ' {2 - 2b) 

with 

D(k) = ai(l-ai)Q 2 + k 2 . (2.27) 

In eq. (2.26) the first factor (with subscript '«') comes from the fermion loop, the second factor 
from gluon production below the fermion loop. For the latter we have made use of the i^ 2 ^4 gluon 
vertex [16, |l~7| ]: 



k 2 + l k 2 \ / k 2 + l' 



(k 2 + l) 2 k 2 ;A(k 2 + i') 2 k 2 /; 



(2.28) 



(this includes the two vertical gluon propagators with momenta k 2 ). A similar expression holds for 
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Next we turn to the four diagrams in Fig. ^ where four gluons are attached to the fermion 
loop. We write : 

/j2i r fly 

^ (XJP ' J ^ {X ^ 1/2 )( T « T ^) (61) ( 2 - 29 ) 
The first diagram in Fig. |||& leads to: 

, {bl) = f l + k!+k 2 1-ki _ kj+k 2 | ki \ / k 2 +l _ k2\ 

1 11 3° ~ U(l + k! + k 2 ) ^DCl-kx) D(k 1+ k 2 ) + ,D(k 1 )JA(k2 + l) 2 k|^ 



r + kl+k2 + (T^-O) (-30) 



ki 


\ / k 2 + l 


D(k x ) 


/A(k 2 + 1) 2 


ki ^ 


\ / k 2 + r 



^(l' + ki+ka) D(l'-ki) D(k 1 +k 2 ) DQl x )J V(k 2 + r) 2 k 2 
Similarly, we find for the other three contributions: 
(T rpi s{b2) = ( 1 - ki - k 2 1 + ki _ ki +k 2 ki \ / k 2 - 1 



VD(l-ki-k 2 ) ' D(l + ki) ,D(k 1 + k 2 ) ' DCki)^ V(k 2 -1) 2 k 2 

( i' + k! + k 2 r-iq k 1 + k 2 , ki \ / k 2 + v 

+ 7T7T7 iTT TTTTT TTT + 7T7T7T ) I 7T HA? ^2 I A^-U 



V^l' + ki + ka) D(l'-ki) ^(ki+ka) DQlJJ , V(k 2 + l') 2 ki 



v 2 



, (63) = _ / l + kx + k 2 1-kj _ kx + k 2 ki \ / k 2 + l _ k, 

1 11 ^ UKl + ki + k^Dfl-ki) ,D(k 1 + k 2 ) + J D(k 1 )/A(k2 + l) 2 k 2 

/ l'-ki-k 2 l' + ki k!+k 2 ki \ / k 2 -l' 



and 



V J D(F-k 1 -k 2 ) + ,D(l + k 1 ) D(k 1 + k 2 ) + D(k 1 )jA(k2-l / ) 2 k 2 ' - (2 - {L>) 



, (M) = / 1 - ki - k 2 1 + ki k x +k 2 kj \ / k 2 - 1 _ k^ 

1 a 3° " U(l-k x -k 2 ) + J D(l + k 1 ) D(k 1+ k 2 ) + ,D(k 1 )JA(k2-l) 2 k 2 

f - ki - k 2 r + kx _ k x + k 2 kx \ / k 2 -r _ k,\ 

^(l'-kx-k^ + ^F + kx) J D(k 1 +k 2 ) + D(k 1 )JA(k2-l / ) 2 HJi" 

The third group in Fig. ||]c consists of diagrams with three gluons being attached to the quark 
loop. We obtain: 

(T . mi s{d) _ 9 ( ki + k 2 _ kx \ / k 2 + l _ koA 
{la ^ l) 2 U(kx+k 2 ) J D(kx)JA(k 2 + l) 2 



l' + kx+k 2 , l'-kx kx + k 2 , ki \ / k 2 + l'_ 

2 \ 



£>(!' + k x + k 2 ) + D(Y-k 1 ) D(k 1 + k 2 ) + Z)(k 1 )JA(k 2 + l') 2 k 2 ' (2 -" >1) 



( ki + k 2 ki \ / k 2 + 1 



\D(k 1 + k 2 ) D(k 1 )J l V(k 2 + 1) 2 k 2 ; ; 

v - ki - k 2 r + kx _ kx + k 2 i kx \ / k 2 -r _ k2 

D(l'-kx-k 2 ) + D(V + k 1 ) D(ki+k 2 ) + J D(k 1 )JA( k 2-l') 2 k 2 



(T T> ^( c3 ) - 9 I 1 + k i + k 2 , 1- ki _ ki + k 2 ki \ f k 2 + l 



D(l + ki+k 2 ) D(l-ki) £>(ki+k 2 ) D(k 1 )J i V(k 2 + 1) 2 k 2 

/ ki + k 2 _ ki \ / k 2 + r _ k 2 



V^(k 1 +k 2 ) D(ki)y , V(k 2 + F) 2 k? 



(2.36) 
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, (c4) _ _ / l-k!-k 2 Hhki _ k 1+ k 2 ki \ / k 2 -l _ k,\ 

1 11 3l) Ua-kx-kaj^^i+kx) D(k 1 +k 2 ) + D(k 1 );^(k 2 -i)2 

ki +k 2 ki \ / k 2 + r k 2 \ 
D(k 1 + k 2 ) D(ki)JA(k2 + l / ) 2 HJi 



(2.37) 



Combining all contributions (a) - (c) we arrive at our final result. The sum of the products 
TuT^ can be written as the square of two identical (up to the primed /-integral) expressions: 

= / l + ki+k 2 ki +k 2 ki -1 _ ki \ / l + k 2 _ ksX 

*' VJ D(l + k 1+ k 2 ) + ^(kx+k,) D(ki-l) D(k 1 )) t v (l + k 2 )2 k 2 J l 

+ (1 -► -1) . (2.38) 

In order to obtain M/, we simply drop the numerators in the first factor (with subscript '«') in 
eq. (2.38): 

/ 1 1 1 1 \ / 1 + k 



D(l + k! + k 2 ) ' £(ki + k 2 ) D(ki-l) D(ki) A(l + k 2 ) 2 \4J t 

+ (l->-l) (2.39) 

which defines 7), and M; then follows from eq. (2.24). Inserting these results into eq. (2.12) - 
eq. (2.15), we arrive at our final expression for the ep cross section. 



3 Impact Parameter Representation 



It is very instructive to analyze the structure of the matrix element Mu or M; in impact parameter 
space. As will be shown below the simplification in impact parameter space arises in form of two 
wave functions of which the first describes the dissociation of the photon into the qq - pair and 
the second the subsequent radiation of the gluon from the quark or the antiquark. The couplings 
of the t - channel gluons generate phase - factors which after some regrouping can be turned into 
effective dipole cross sections. For the numerical analysis later on we have, of course, to use the 
momentum representation. 

We start by taking the Fourier transformation of Tu : 



J d 2 ki d 2 k 2 e ikl ' ri+ik2 - r2 T a 



(3.1) 



Vff(ri) 



2 k e ik 2 -r 2 



| e -i(k 2 +l) ri _|_ e -ik 2 -ri _ g-il-ri _ ]_ j | 



+(!->-!)] 

W (r 2 - ri) {e^ 2 + e'^^ - e -* ri - l} + (1 
- <^(r 2 ) {e"^ r2 + e -^ r2 - ri ) - e iZ ' ri - l} + (1 -» 



k 2 + l 

(k 2 + l) 2 

-1) 



k 2 , 

K 2 > I 



where \E ,T and ip are defined as 



*f(r) 



d 2 k e ikr 



/ 



d 2 k e lk r h 



D(k) 

k 2 



|r| 
r/ 



-2m -± \la{l - a)Q 2 Ki(^/a(l - a)Q 2 r 2 ^ 
-2iti 



(3.2) 
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If we now introduce the effective dipole cross-section as 

d 2 l 



(TIT) 



J |i F(x F , I 2 ) (l - e^ r ) (l - e-^ r ) (3.3) 



we finally obtain for Mu, the Fourier transform of Mu'. 

Mu = *f(ri) Wi(r 2 )-i Pl (r 2 -r 1 )} l [a(r 2 ) + <r(r 2 - n) - a( ri )] . (3.4) 

The vector ri denotes the separation of the quark and antiquark, r 2 the separation of the quark 
and gluon and r 2 — ri the separation of the antiquark and gluon. The result in eq. (|3.4j) is the same 



as in ref. [18|. One derives the corresponding expression for Mi by substituting the transverse wave 
function *f? T by the longitudinal wave function ty L , 

* L (r) = J d " k e*" -^y = 2vr K (y/a(l - a)Q 2 r 2 ) . (3.5) 

The impact parameter representation suggests the intuitive interpretation that following the initial 
dissociation of the photon, which is described by ty T)L , each of the quarks radiates off a gluon 
represented by (p. After the creation of the qqg-st&te the interaction with the target processes 
pairwise, i.e. each of the pairs qq, qg and qg gives a separate contribution as indicated by the 
different arguments of the effective dipole cross-section a. As was pointed out in ref. flql , however, 
the color structure is not quite consistent with the previous interpretation. In particular the 
interaction of the qq as part of the qqg-Hnal state is color suppressed by powers of N c . The reason 
for having a non-suppressed contribution from the qq - pair is hidden in the requirement to include 



the interaction of the 'renormalized' qq - pair [18]. In other words, the expression in eq. ( |3.4|) also 



includes a contribution where the interaction with the qq takes place before a gluon is emitted. 



This contribution is not color suppressed. The result eq. (3^1) has been derived from Feynman 
diagrams which automatically take into account all necessary configurations. 
Taking the square of the amplitude in impact parameter space one finds 

r 2 

M U M U = (2ir) 2 ^(rx)] 2 ^ 1 [a(r 2 ) + a(r 2 - n) - a(n)] 2 . (3.6) 



l'2-ri 



The expression r s p i s characteristic for the dipole formalism |T{J. The exact equivalence of 

eq. ( |3.6| ) and the corresponding double - dipole scattering still needs to be established. 

Assuming strong ordering in the separation of the quark - antiquark - pair and the separation 
of the quarks and the gluon, |ri| <C |r 2 |, one obtains the simple factorized form 

MaMa = (2vr) 2 r 2 |^( ri )| 2 4^ . (3.7) 



4 A Few Special Kinematic Regions 

In this section we consider a few kinematical regions of interest for which analytic expressions can 
be obtained. We will average over azimuthal angles. As a result, we need to consider only eq. (2.12) 
and eq. (2.14). Eq. (2.15) would require a discussion similar to eq.(2.14). 
First we study the 'hard' region: 

Ql < l 2 < k 2 ~ D(ki) ~ D(ki + k 2 ) ~ k 2 (4.1) 
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where Qq denotes a hadronic scale which separates the hard and the soft region. For simplicity 
we assume that all relevant scales in the problem are of the same order of magnitude, k. In the 
7* — JP CM-system the gluon jet is close to the direction of the Pomeron, and the two quarks 
form a broader jet in the opposite direction. For Tu in eq. (2.38) we find (after averaging over the 
azimuthal angle of 1): 



Tu 



1_ 2 

4 



,(ki + k 2 )i(ki + k 2 ) 



k 2 , V£>(ki + k 2 ) 



+ 



,(ki);(k 



1 im 



£)2( kl + k2 ) 1 D(ki) 
Using eq. (2.24) and eq. (2.25) we obtain for MuM' u : 

MaM', = 2 (/ ' ' 



^ 2 (ki 



'ml 



(k 2 ) m (k 2 ); 



k| 



(4.2) 




L»2( kl ) 

[kl • (kl + k 2 )] 2 - 



(4.3) 



k 2 (k 1+ k 2 ) 2 



\D 2 {k! +k 2 



,D 2 (k 1 +k 2 ) J D 2 (k 1 ) 



X]pg(xip,k 2 



The transverse cross section is obtained by inserting eq. (4.3) into eq. (2.12). The most remarkable 
feature of this result is the appearance of the square of the gluon structure function: because of the 
prefactor l/(k 2 ) 2 , one expects the 'hard' region of large k 2 to be suppressed. But as a result of its 
rise at small xp, the gluon structure function provides an enhancement factor of the hard region 
which increases with decreasing xp. 

As a special case of the hard region we consider the case of strong ordering of the transverse 
momenta, k 2 , <C k 2 . Then eq. (4.3) simplifies into 



Tu 



l 2 



,(kiMk 



1 Im 



k 2 \D(k!) 



D 2 (k 1 ) 



5ml 



(k 2 ) m (k 



2 1 



k 2 

K 2 



(4.4) 



With the approximation of eq. (2.19) we obtain for MuM-,: 



MuM^ = 8 



Q 4 + m 4 m 2 \ 1 
(Q 2 + m 2 ) 4 I kf J (kfF 



xpgixjp^l) 



(4.5) 



It might be interesting to note that an early approach [20] to diffractive deep inelastic scattering 
which not only considers the radiation of one but many gluons was based on the special case 



discussed above. An extension of these results towards large (3 can be found in ref. [21] which in 
our notation reads 



MuM' u 



Q 4 + m 4 
\Q 2 + m 2 ) 4 



m 



M 2 



m 



, 4 /M 2 + 2m 2 + 3Q 2 \ 2 



1 



M 2 + Q 2 



M 2 + Q 2 



(k|) 2 



xpgixp^l) 



The main observation here is the strong suppression of MuM^ when m 2 becomes of the order of 
M 2 . This fact is also implemented in the parameterization of Q. 

Next we say a few more words about the small-k| region. As we have noted already after 
eq. (4.3), the region of small k 2 is expected to give an important contribution to the integrated 
diffractive cross section. On the other hand, if k 2 becomes smaller than, say, the hadronic scale 
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Qq, perturbation becomes unreliable and we cannot apply our formulae. For practical purposes, 
however, it would be attractive to find an extrapolation of our jet cross section into the region 
where the gluon jets becomes soft. Following the discussion of the qq final state in diffraction 
dissociation in || we shall present an extrapolation of our perturbative calculation, which at small 
k| leads to the Ingelman-Schlein picture of the Pomeron structure function. First, at small k|, the 
approximation l 2 -C k| will no longer be valid; l 2 can become of the same order as or even larger 
than k 2 .. For simplicity we do not consider the most general case but still retain the approximation 
l 2 , k| <C k 2 . After integrating over the azimuthal angle of 1 we find for Tu : 

Tu = 2 {DM ~ 2 D 2 ^) ) V ml ~ 2 —£ ! H 6(1 -k 2 ) + ^e(k 2 -l)J (4.6) 

which is consistent with the result found earlier in ref. pi]]. Before inserting eq. (4.6) into eq. (2.29), 
we have to discuss the two gluon amplitude J-(xjp,l 2 ). As long as the momentum scale I 2 at the 
upper end is large, it can (to a good approximation) be identified with the unintegrated gluon 
structure function (see eq. (2.25)). Now we will try to find an extrapolation into the soft k 2 region 
which, at low k 2 , smoothly turns into the Pomeron structure function picture. An easy way has 
been sketched in ||, and here we outline the argument in somewhat more detail. We make the 
following simple double Mellin transform ansatz for J-: 

^,l 2 ) = ^/^/^(^) (4.7) 

where the integration contours run along the imaginary axis (— 1 < Refi < 0, x(/i) < Reco), and 
())q denotes the overall normalization constant. If we would perform the [i integral first by picking 
up a pole at /i = \~ 1 (uj) = 7(0;), we obtain the anomalous dimension. However, here it is more 
convenient to model the function x(p) an d to do first the u integral. Let us assume that x(^) is 
positive and becomes large both at = —1 and /U = (similar to the BFKL characteristic function). 
We then obtain: 

Tixjp, l 2 ) = o i / 7T-- exp ( x0«) log—+fi log ] (4.8) 
J 2m \ xjp / 




Combining this ansatz with eq. (4.6) we write the I integral as: 

" l!Sr ( l ' 2 Y [ ( 1/2 - 1) + ^©(1 - 1*)] (4-9) 

Inserting the k^-dependent prefactor of eq. (4.9) into the integral of eq. (4.7) and performing a 
stationary phase analysis, we find that for Qq <C k 2 , the stationary point moves towards fj, = — 1, 
whereas for the opposite case k 2 , <C Qq the stationary point is near fi = 0. In order to evaluate the 
l' 2 integral in eq. (4.8) near \x = — 1 or p = we have to specify the lower limit of integration. In 
the first case (large k|, /i — > —1), the first term in eq. (4.9) stays constant whereas the second part 
gives a logk^/Qg from the region Qq/^I < l' 2 < 1- In the second case (small k|, /U — > 0) the first 
term gives a logarithmic enhancement logk 2 /(5g whereas the second term stays constant (which is 
equal to 1 if the I' 2 integral extends down to 0). As a result we obtain for the l 2 integral: 

,2 1,2 



d>(kf,k^,xjp) 




(4.10) 
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Together with eq. (4.6), we obtain for MuM^: 

. 8m 4 Q 4 + m 4 

Inserting this into (2.13), and making use of (2.20) we find the expected result: the transverse cross 
section belongs to leading twist, and the integral over k 2 diverges logarithmically. 

A nice feature of this extrapolation into the soft non-perturbative region is the following. As 
we have described above, the two cases (large and small k 2 ) correspond to different regions of the 
l 2 integral: the first case is the usual strong ordering with Qq < l 2 < k|. In the second case l 2 
prefers to be larger than k 2 .: in this region we simply drop k 2 , in comparison with 1 and ki. As a 
result, the non-planar diagrams in Fig. |] become small in comparison with the planar ones. What 
is left can be interpreted as a 'gluon component of the Pomeron structure function'. This has to be 
compared with the first case, where all diagrams are equally important, i.e. the Pomeron interacts 
with the whole diffractive system. 

Next we consider the case where the gluon stays hard (large k2) and one of the quarks becomes 
soft, i.e. ki ~ or ki + k2 ~ 0. Let us assume ki to be small (k 2 < l 2 < Qq) and go back to 
eq. (2.38). With l 2 > k 2 only the third term in the first bracket contributes and we find after 
integration over the azimuth angle: 



,(k 2 Mk 2 )A J_ 2 
k 2 J D(l) 



T il =(5 il -2 " '>l\ ^ )-^ . (4.12) 



Using our ansatz for the unintegrated structure function eq. (4.7) with the soft extrapolation /x — > 
we can perform the 1-integration 



$(/ci,xjp) 



m d\ 2 ( i 2 V i 2 



k? l 2 \QIJ DO) 

( D(Qo) 
U(k0 



... m i (-3) 



and finally arrive at 

MuMit = 1> 2 ^ ■ ( 4 - 14 ) 

It is important to note that this kinematical limit yields a contribution which breaks (collinear) 
factorization. The result in eq. (4.14) is similar to the 'super-hard component of the Pomeron' 
introduced in ref. [14]. In a scenario where the photon virtuality Q 2 is smaller than the gluon 



momentum k|, one would expect from the quark loop a large logarithm of the form log^l/Q 2 
if factorization worked. The potential softness of one of the quarks which controls the scale l 2 in 
eq. (4.13) is responsible for the lack of such a hard logarithm. In other words the soft part of 
the process decouples from the jets in the final state. With regard to the inclusive diffractive cross 
section, the hard gluon jet combined with a soft quark gives only a rather small fraction of the total 
contribution. The leading contribution comes from the scenario discussed earlier where the gluon 
becomes soft and the quarks are hard. In this case we have the familiar strong ordering situation 
which provides a hard logarithm log(Q 2 /k 2 ) in accordance with factorization for diffractive deep 
inelastic scattering [15] and the notion of a Pomeron structure function. Factorization is still 
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violated by contributions with k| bigger than the assumed factorization scale. Since the mean 
value of k 2 , is close to Qq, as can be inferred from eq. (4.10) and eq. (4.11), one can argue that for 
any scale well above Qq, factorization should hold. 

Finally we return to our 'hardness condition' in eq. (2.2). For the gluon jet, its transverse 
momentum square k 2 , coincides with the virtuality of the vertical gluon line in Fig. ||a: for small 
k 2 . this 'Pomeron structure function' diagram becomes more important than the other non-planar 
diagrams; the boundary between hard and soft physics is defined by the off-shellness of the gluon 
emitted from the qq pair. Correspondingly, for the quark-antiquark pair it is the virtuality of the 
quark (or antiquark) before it emits the gluon. With the approximations discussed in the beginning 
of this section, the requirement that the four momentum square of this quark has to be large then 
leads to the last two equations in eq. (2.2). 



5 A Numerical Study 

In this section we study the specific signature of the process e + p — >q + q + g + p' with the cross 
section calculated in the previous sections. We work in the j - IP center of mass system, and the 
gluon will be restricted to the forward hemisphere of the proton (Pomeron). We have implemented 
the formulae of eq. (2.12) - eq. (2.17) into a full hadron level Monte Carlo program | |22| |. To 
avoid soft divergences of the matrix element we impose a cut in \t\ = \(q — Pq) 2 \ > 2 GeV 2 and 
H = I (q ~ Pq) 2 \ > 2 GeV 2 , where q, p q = q + r — k\ — A;2, and Pg = k\ are the four momenta of the 
photon, the upper outgoing quark and the lower outgoing antiquark (Fig. |||), respectively, and a 
cut on the transverse momentum squared of the gluon k 2 , > 2 GeV 2 . The coupling a s is kept fixed 



at a s = 0.25. For the gluon structure function we use the NLO parameterization of GRV [23|. We 
divide this numerical study into two parts. In order to analyze the general features of the three- 
parton final state we first work in an "ideal" (i.e. truly asymptotic) kinematic environment where 
the invariant mass M is much bigger than all transverse momenta (Q 2 = 100 GeV 2 , M 2 = 900 
GeV 2 , (i.e. = 0.1), xjp < 10~ 2 ). In the second part we turn to HERA kinematics. 



5.1 Asymptotic region 

First we have to define the region of validity of our calculation. In deriving our cross section formulae 
we have made several approximations. Regge-kinematics requires that «2 "C <*!■> 1 — a\. In our 
numerical analysis we therefore impose a kinematical cut and require that 02 < g«i> |(1 — 
Second, we expect that for large M 2 the mass of the qq system will be small: in order to enforce 
this condition we further demand that m 2 g < \M 2 . 

We begin with pt spectra (Fig. ||). A fit to a power-like behavior (l/p^) p (for p\ < 10 GeV 2 ) gives 
p = 0.82 and p = 1.63 for the quark and gluon transverse momenta, respectively. This behavior 
is not far from the naive estimate based upon eq. (4.5): together with the phase space factors in 
eq. (2.12) - eq. (2.13) one expects p = 1 for the quark momentum k 2 , and p slightly less than 2 for 
the gluon momentum k 2 .. A very striking feature is the fact that for a rather large fraction of events 
the ordering condition k 2 < k 2 , (ki + k2) 2 is not fulfilled: we find that only 60% of the events 
satisfy this condition. This demonstrates that the simple boson-gluon picture in which the gluon 
is emitted with a transverse momentum much smaller than that of the two quarks may be very 
unreliable. On the other hand, for almost all events the square of the gluon transverse momentum 
is less than the virtualities of the exchanged quarks (k 2 . < \t\, \u\). 



16 





Figure 7: The distribution of thrust T. 



Next we look into the spatial distribution of the three-parton final state. From the experimental 
point of view it is convenient to begin with the thrust distribution: the plot in Fig. [?] shows that 
a large fraction of the events has a two-jet like structure (T > 0.9). However, there remains also 
a sizable fraction (about 20%) of three-jet events for which T < 0.9. To analyze these samples 
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Figure 8: The distribution of the opening angle A of the qq pair for values of thrust 
T > 0.9 (two jet configuration). 

in more detail, we begin with the two-jet events and we show in Fig. |8| the distribution of the 
opening angle A of the qq pair: the distribution peaks at both small and large opening angles. We 
divide the events into two sets: events with large opening angle (A > tt/2) and those with small 
opening angle (A < 7r/2). In the first case we find (Fig. ^) a large asymmetry in the longitudinal 
momenta and a big invariant mass of the two quarks. The quark which moves with maximum 
momentum (~ 15 GeV) in the photon direction forms one jet whereas the second quark moves 
into the opposite direction. The second quark and the gluon form the second jet sharing the jet - 
momentum. In Fig. [l0| a the correlation of the longitudinal quark and gluon momentum for T > 0.9 
and A > tt/2 (large opening angle) is shown. The upper band contains events where the quark 
moves isolated in the photon direction whereas the antiquark moves along with the gluon in the 
Pomeron direction. The lower band contains events where the quark and the gluon form a jet 
sharing the jet-momentum. 

For those events with small opening angle, on the other hand, the momenta of the two quarks 
are much more symmetric. In this case, the two quarks form a jet opposite to the gluon (Fig. |9p. 

In Fig. |i~0| fr we show the quark and antiquark longitudinal momentum for T > 0.9 and A < tt/2 
(small opening angle). Both momenta are strongly correlated. 

In all these cases the jet axis (i.e. the direction of the gluon or quark) lies mainly in the Pomeron 
direction ("aligned gluon configuration"): Fig. [l2| (full curve) shows the distributions of the two 
- jet events. Compared to the full set of events (Fig. [l^ dotted curve) one notices a slightly steeper 
decrease. 

For the three jet events (T < 0.9) the opening angle A ranges between 1.2 and 2.2 (Fig. ^), i.e. 
we have a star - like configuration. Again, the gluon prefers to be in the forward direction (Fig. [l2|, 
dashed curve). However, compared to the two - jet events, the fall - off at large k| is less steep. In 
other words, events with a large gluon transverse momentum are more likely to be three - jet like 
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Figure 9: The difference of the longitudinal momenta of the two quarks along the z 
direction (A(p z )) as a function of the opening angle A of the qq pair for thrust T > 0.9. 
The second plot shows the correlation of the gg-invariant mass and A. 

than those with a small transverse momentum. A few numbers: for k| < 10 GeV 2 the ratio of 3 - 
jet to 2 - jet events is about 0.2, whereas for k 2 . > 10 GeV 2 this ration changes to 0.47. 

5.2 HERA Region 

After this more general study we now turn to the HERA region. The results of the previous study 
should apply as long as the diffractive mass is sufficiently large: only under these conditions we 
have a sufficiently large range in transverse momenta. In HERA kinematics this requirement forces 
us into the (difficult) region of very small (3: for Q 2 = 5 GeV 2 and M 2 = 500 GeV 2 , we would 
need (3 = 0.01. In a first step we have repeated our above analysis for this region, and we have 
found that, indeed, the previous conclusions on our model still apply. In this region, however, 
we expect secondary exchanges to become important, and the theoretical interpretation of a jet 
analysis becomes more complicated. 

A more realistic /3-region is j3 = 0.1. In the HERA region of not too large Q 2 we are then 
limited in M 2 . This restriction, together with the lower cutoffs on the transverse momenta and 
the virtualities of the partons, lead to a severe limitation of the phase space for the transverse 
momenta of the jets. In our numerical analysis we have chosen Q 2 = 5 GeV 2 and M 2 = 45 GeV 2 . 
We do not impose a cut on m 2 g or a\ and Q2- Such a cut would suppress the cross section by a 
factor of around 10 rendering the analysis meaningless. Our strategy will be a comparison with the 
asymptotic situation. From this we can deduce how far we can trust our asymptotic formulae for 
the HERA kinematics. 

We again begin with the fcj-spectra (Fig. [14|) for the proposed HERA kinematics, and we fit 
the power p. For the quark-spectrum we find p = 1.8 and for the gluon-spectrum p = 2.3. Both 
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Figure 10: a. The correlation of the longitudinal quark and gluon momentum for T > 0.9 
and A > tt/2 (large opening angle), b. The longitudinal momenta of the quark and 
antiquark for T > 0.9 and A < tt/2 (small opening angle). Both momenta are strongly 
correlated. 



distributions have become steeper as compared to the asymptotic case which is mainly due to the 
stronger restriction in the phase space. The available energy is approximately a factor 10 smaller 
now. We also learn that a measurement of the /cj-spectrum cannot directly be related to the 
analytic estimates we have done in the previous section. In order to test to what extent the usual 
assumption of strong ordering between the transverse momenta of quarks and gluons is satisfied we 
compute the fraction of the cross section where k 2 . < min(|i|, \u\). we find that approximately only 
2/3 of the cross section satisfies this constraint, whereas 1/3 lies outside this region. Next we have 
a look at the thrust distribution (Fig. 15). It is much broader than in the asymptotic case. Since 
M 2 is no longer that large, the relative proportion of transverse to longitudinal momenta is much 
higher than before, and therefore leads to an increase of three-jet events (around 60% of the events 
have a thrust smaller than 0.9). If we ask for two-jet events, i.e. thrust T > 0.9, then we find again 
two configurations (Fig. [jl]), one with rather small opening angle between the quarks and a second 
with large opening angle. The distribution of Ap z versus opening angle A looks similar to the 
corresponding plot for asymptotic energies which makes us believe that the asymptotic formulae 
also work properly for HERA kinematics. 

Table 1 shows the integrated qqg diffractive jet cross sections. For comparison, we also present 
dijet cross sections from qq production. We have integrated over the kinematic region 10 GeV 2 
< Q 2 , 50 < W < 220 GeV, x F < 10~ 2 , where we have integrated over t, assuming an exponential 
t dependence. In the first row we present cross sections for the case where all parton k 2 are larger 
than 2 GeV 2 in addition to > 2 GeV 2 , and in the second row the corresponding cut is 

5 GeV 2 . For comparison, in the same kinematic region the inclusive diffractive cross section (for 
three flavors) is estimated to be ~ 2.5 nb, i.e. for the lower kt cut the combined qq and qqg jet rate 
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Figure 11: The distribution of the mass of the pair for thrust T > 0.9 

amounts to approximately 8% of the diffractive cross section. A look at the kt spectra indicates 
that the qqg jet rate strongly increases if we lower the cutoffs (entering a region where our formula 




Figure 12: The transverse momenta of the gluon for events with T > 0.9 (solid line), for 
T < 0.9 (dashed line) and for all events (dotted line). 
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Figure 13: The difference of the momenta of the two quarks along the z direction (A(p z )) 
as a function of the opening angle A of the qq pair for thrust T < 0.9. 

requires a modification of the gluon structure function) : this indicates that the largest part of the 
qqg cross section might come from an intermediate region where, in particular, the gluon kt is larger 




Figure 14: The distribution of transverse momenta for quarks and gluons. 
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Figure 15: The distribution of thrust T at HERA energies. 
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Figure 16: The opening angle A for T > 0.9 at HERA energies. 



than the soft Pomeron scale but lower than the cutoff values used in our analysis. 

As in the case of diffr active production of two jets we expect the qqg jet cross section to 
rise in 1/xjp; for the special case k| <C k 2 such a behavior follows immediately from eq. (4.5): 
a ~ [xipg(x]p, /U 2 )] 2 . In Fig. [l?] we show the cross section as a function of xjp at fixed (3. In the 
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kinematic cut 


dijet 


qqg 


2 GeV 2 
5 GeV 2 


117 pb 
55 pb 


66 pb 
9.2 pb 



Table 1: Integrated cross sections for diffractive qq and qqg production in the kinematic 
region denned by: 10 GeV 2 < Q 2 , 50 < W < 220 GeV, x F < 10" 2 . The t dependence is 
assumed to be exponential and the t integral is performed. The kinematic cut requires all 
parton k 2 and in the case of qqg also \t\, \u\ to be greater than the value specified. 

-2 

x 10 




pom 



Figure 17: The cross section as a function of xjp at fixed (3 = 0.1. The cross section is 
integrated over t. 

region /3 = 0.1 (i.e. large diffractive mass) we expect our formula to be applicable. One can see 
the steep rise of the cross section at small xjp, a similar behavior as observed already in the dijet 
case. This rise is in qualitative agreement with the approximation eq. (4.5). As we have indicated 
in the introduction, we expect the qqg diffractive final state to dominate in the region of large 
diffractive masses (small j3). In Fig. |l^ we show, for different j3 bins, the integrated cross sections 
of the two diffractive processes qq and qqg production. One sees that the qqg becomes dominant 
approximately for (3 < 0.2. It is interesting to note that this (3- value is consistent with the fit to the 
diffractive cross section in ||: this fit also suggests that the qqg contribution begins to dominate 
for (3 ~ 0.2. In Q, the solution with a large 7- value is closer to our estimate than the small-7 
solution. 

Next we study in more detail the specific properties of the qqg system, in particular the angular 
distribution of the three jets. We begin with a scatter plot of the angles of the quarks, keeping the 
angle of the gluon fixed. Quarks angles are defined as angles between quark and photon directions. 
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Figure 18: Diffractive jet cross sections da/d(5 at t = for 0.01 < y < 0.9, Q 2 > 5, and 
xjp < 0.01. All fcy's are restricted to be larger than 2 GeV 2 . The solid lines shows the qqg 
final state, the dotted one the final state qq. 

In Fig. |19] we show the results for three different regions of the angle 62 of the gluon. According 
to eq. (2.22), this angle is related to the transverse momentum k2 of the gluon. At small gluon 
angle (direction of the gluon close to the proton (Pomeron) direction, Fig. |i~9|a,), in most of the 
events both quark and antiquark are moving not far from the photon direction (i.e. both their 
cost's are close to 1). However, there are also events in which one of the quarks stays close to the 
photon direction, whereas the other one sticks out at a larger angle. At large angle of the gluon 
(gluon direction orthogonal to the photon-Pomeron beam axis, Fig. |T~9]c. ) , on the other hand, the 
quarks tend to move in opposite directions, one in the photon direction and the other in the proton 
direction. Correspondingly, at small gluon angle the opening angle A between the quarks prefers 
to be around 7r/3 (Fig. |20"|a.), whereas for large gluon angle (Fig. |20|c.) A peaks at a value slightly 
above 2ir/3. In all these discussions it should be made clear, that case (c.) lies at the edge of the 
region of validity of our approximation. A reliable calculation of this kinematic region requires a 
complete NLO calculation of diffractive qq and qqg production, in particular virtual corrections to 
qq production and real gluon emission from quark lines. 

From these considerations we arrive at the following picture: 

• At small transverse momenta of the gluon, we have mainly two-jet configurations where the 
quarks travel in approximately the same direction with an opening around 7r/3, carrying 
similar longitudinal momenta. But the two quarks can also appear in a more asymmetric 
configuration, where one of the quarks has a much larger longitudinal momentum compared 
to the other quark, and moves in approximately the direction opposite to the gluon. The 
quark with the smaller longitudinal momentum then sticks out from the z-axis. 

• In the other extreme, where the gluon is emitted perpendicular to the z-axis, we observe a 
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Figure 19: The 9 angle of the quark versus the 9 angle of the antiquark. (a.) for —1 < 
cos(9 gluon ) < -0.8, (6.) for -0.6 < cos(6 glu(m ) < -0.4, (a) for -0.2 < cos(9 gluon ) < 0.0, 
at = 0.1 and Q 2 = 50 GeV 2 (M\ = 45 GeV 2 ). 

'Mercedes — star" like configuration, in which both quarks carry longitudinal momenta of 
similar value, but in opposite direction, and the cosine of the opening angle lies between two 
and three. 

• At medium transverse momentum k2 of the gluon, all configurations discussed above, are 
possible. 

6 Conclusions 

In this paper we have presented an analysis of hard diffractive qqg jets with large transverse mo- 
menta. Our study is motivated by the observation that in DIS diffractive dissociation, the Pomeron 
intercept lies above the value seen in purely hadronic interactions: this suggests that the diffractive 
final state contains a rather large 'hard' component. A natural candidate for hard final states are 
jets. The investigation of diffractive qqg production presents a generalization of diffractive dijet 
production ||. Our calculation is restricted to the low-/? region (large diffractive masses), and 
we have been working in the leading-log 1/(3 approximation. At the same time, the transverse 
momenta of the outgoing partons are not restricted by a strong ordering requirement. In the low- (3 
region the emitted gluon carries a small momentum fraction of the incoming photon, i.e. in rapidity 
it is closer to the proton than the two quarks (Fig. |3|). 

The analytic expression for the cross section formula in impact parameter space is remarkably 
simple. It illustrates that the concept of a photon wave function holds even beyond the approxi- 
mation, where the transverse momentum of the emitted gluon is much softer than the quarks. For 
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Figure 20: The opening angle A between the quarks, (a.) for — 1 < cos(6 g i uon ) < —0.8, 
(6.) for -0.6 < cos(9gi uon ) < -0.4, (a) for -0.2 < cos(0 gluon ) < 0.0, at (3 = 0.1 and 
Q 2 = 50 GeV 2 {M\ = 45 GeV 2 ). 

special kinematic limits we have obtained even simpler expression for the cross section formulae. 

In a first exploratory numerical study we have calculated jet cross sections. Depending upon 
the lower limit of the transverse momenta, the combined qq and qqg jet cross section can be as 
large as 8% of the inclusive diffractive cross section. We have also looked into the spatial distri- 
bution of the three-parton final state. If (in the 7* — IP CM-system) the gluon jet is close to the 
Pomeron direction, the quark - antiquark pair mainly forms a single (although somewhat broader) 
jet opposite to the gluon jet. In the extreme case where the emitted gluon is perpendicular to the 
photon - Pomeron axis, a clean three-jet configuration emerges. However, in this region the cross 
section will be small. For a diffractive jet analysis, it is important to keep in mind that quite a 
large fraction of the qqg final states will appear as a two-jet configuration, and it will be difficult 
to separate them from the qq dijet final states. 

Experimentally the analysis of two or three jet events is not easy. A more promising way may 
be a thrust or Et analysis. To this end, it would be very desirable to have a complete NLO 
calculation of diffractive final states with two or three jets. 
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